m 
O 

o 



On Gauged Renormalisation Group 
Transformations of Lattice Fermions 



o 
o 

(N 

5-1 . 

Ph. 

<^ Artan B 01191 

oo 



University of Tirana 
Department of Physics, Faculty of Natural Sciences 
King Zog I Boulevard, Tirana, Albania 
borici @fshn. edu. al 



Abstract 

We construct a hierarchy of lattice fermions, where the coarser lattice Dirac operator 
is the Schur complement of the block UL decomposition of the finer lattice operator. We 
show that the construction is an exact gauged renormalisation group transformation of the 
lattice action. In particular, using such a transformation and the QCDLAB tool, it is shown 
how to implement the Ginsparg-Wilson strategy for chiral fermions in the presence of a 
dynamical gauge field. The scheme allows, for the first time, a full multigrid algorithm for 
lattice quarks. 
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1. A renormalisation group transformation for quadratic actions, as it is the case of lattice 
fermions, is a simple Gaussian integration, 



det D e -^ bSbh1pb 
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where D and Sbb are Dirac operators on the fine and coarse lattices, D bb is a naive Dirac operator 
on the coarse lattice, B, B are blocking operators, and, by evaluating the right hand side, one 
can show that 

D = D + BD bb B, S bb = D bb - D bb BD- l BD bb . 

This general transformation can be made concrete by the following prescription: 
1. Partition the Dirac operator as a 2x2 block operator in the form, 



D = 
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2. Define the gauge covariant blocking kernels, 



B =l0 D^D br ), B 
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In this case, one can show that: 



det D = det D bb det D rr , S bb = D bb - D br D r ^D rb . 

On the other hand, a block UL decomposition of D gives 

f > 

hb DbrDrr 

S bb being the Schur complement. Therefore, we may conclude that: 

i) the coarse lattice Dirac operator is the Schur complement S bb of the block UL decompo- 
sition of the fine lattice Dirac operator; 
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ii) the fermion measure of the original theory is det D = det Sbb det D rr ; 

iii) if D is the continuum Dirac operator, the Ginsparg -Wilson relation reads HI : 

{Ts,^ 1 } = {75,^}. 

2. In their paper, Ginsparg and Wilson note that "any h, in particular a fixed point h 
approached after many iterations of the block spin transformation is said to be chirally invariant 
if it satisfies" the relation, u They conclude the paper by noting: "Finding a way to go ahead 
and actually gauge a symmetry present only in remnant form stands as a further challenge" . 

Since then, a lot of efforts have been devoted to solve the Ginsparg-Wilson relation: the 
overlap fermion is an exact solution to this relation 013), whereas domain wall fermions EMU 
and classical fixed point actions (6l lead to approximate solutions. 

Since a fixed point approach in the presence of gauge fields will be a very complicated 
function of the original fine lattice operator, we propose a hybrid approach: trading the Dirac 
kernel of the one step gauged renormalisation group, or possibly its approximation, for the 
kernel of the overlap or domain wall fermions. With one step only, we retain simplicity and 
hope to go away from the Aoki phase. In that case, we would get domain wall fermions with 
smaller residual mass and overlap fermions with better localisation properties. 

Wetestthe idea using the QCDLAB tool and the new functions, Permutat ion_block, 
and Schur.complement. Permutat ionJolock returns a permutation operator P(p), 
which exchanges the rows of D, whereas Schur.complement is used to compute Sbb- The 
permutation p is such that the coarser lattice sites are labelled first, as is shown in Figure [TJ 

We consider two types of permutation: 



'in their notations the relation is {75, h x } = {75, a 1 }. 
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Type II permutation 


6 


7 


8 


5 


6 


7 


8 


9 


10 


11 


12 


1 


2 


3 


4 


1 


2 


3 


4 



Figure 1 : Example of a site permutation on a 4x4 lattice. Coarser lattice sites are labelled by boldface font. 

• Type I permutation: block the lattice in all directions as shown in the upper panel of 
Figure 1. 

• Type II permutation: block the lattice in all directions, one direction at a time, as shown 
in the lower panel of Figure 1. 

Since to every permutation there is an inverse permutation, we have in all four types of permu- 
tations, pi, pi i, pj ,pJi, the corresponding permutation operators being Pj, P n , Pf, Pfj. 

Using these permutations, we have computed four different Schur complements in a SU(3) 
background field using the plaquette action at = 5.4 on a 8 s lattice. At this coupling, the 
spectrum of the Wilson operator should have no spectral gap in the real axes [[8]|. In Figure 
|2] we show the effect of one step renormalisation group transformation on the spectrum of the 
Wilson-Dirac operator. 
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Figure 2: The spectrum of the blocked Wilson-Dirac operator using various permutations of lattice sites as 
described in the titles of the plots. 
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In the case of type II inverse permutation, one can see the emergence of an "owl eyes" plot. 
If such a picture is generic for this particular permutation, we could define the shifted Schur 
complement to be the domain wall or overlap kernel: if c is the centre of the right "eye", the 
kernel of the chiral operator will be cl — S bb . In this case, the benefits are threefold: 

1. better chiral properties/localisation of the domain wall/overlap operator; 

2. faster inversions of the kernel; 

3. better scaling properties of the lattice Dirac operator. 

However, this construction could be computationally expensive: in order to compute the 
Schur complement one has to invert the D„. matrix. In fact, as it will be clear below, a Schur 
complement approximation to the exact one will do. 

3. In principle, one can stay with the exact Schur complement if there is no need to iterate 
the renormalisation group transformation. One exmaple is the domain decomposition approach 
of Liischer He treats exactly the Schur complement and gets excellent results for the two- 
level algorithm. However, if one wants to iterate the scheme to full multigrid one has to rely on 
some approximation of the Schur complement. The reason is twofold: one would like to retain 
the sparsity of the operator and keep the computational overhead under control. 

Earlier papers on Schur complement approximation have ignored the D rr matrix altogether 
IfTOl . Ignoring D rr has the advantage of conserving the sparsity pattern of the coarser operator. 
This is similar to Migdal-Kadanoff approximation, which according to Creutz, "its primary 
drawback lies in the difficulty of assessing the severity of the approximations involved" tfTTI . 
On the other hand, an exact Schur complement may not be practical. The compromise is to 
allow the appearance of one power of D rr . 
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Figure 3: The spectrum of the Schur complement (circles) and its approximation (crosses). 



The Schur complement approximation would be in the spirit of the renormalisation group 
transformation if the coarse lattice operator inherits basic properties of the fine lattice opera- 
tor. In case of the Wilson operator one would like to conserve the positive definiteness, the 
covariance of the hopping matrix, and perhaps 75-Hermiticity. Reusken shows that employ- 
ing point-Gaussian elimination on a weakly diagonally dominant M-matrix, which is similar to 
what we have, one gets a stable Schur complement approximation [fT2l . To the second order, 
his result reads: 



where d rr is the diagonal of D„. and d rr is a diagonal matrix, its entries being the sum of D„ 
rows. Using this approximation, one can compare the properties of S&& and Sbb'- the spectrum 
of Sbb is again an "owl eyes" plot, as in Figure [31 whereas the cost of multiplication by Sj,b is 
approximately the same as the cost of multiplication by D. 

4. In fact, if there is a place to test the Schur complement approximation, this is the full 
multigrid algorithm. Inheritance or stability is of great importance here, as well as regularity, 
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Figure 4: The spectrum of S^Sbb using variuos permutations of lattice sites as described in the titles of the 
plots. 

i.e. how close an approximation is to the exact Schur complement. Without regularity, it would 
be impossible to iterate the two-grid algorithm to a full multigrid. From Figure H] we see that, 
indeed, such an approximation gives a spectrum of S^Sbb clustered around one. Moreover, we 
observe that the type II and type II inverse permutations give excellent results. 

Using such approximations of the Schur complement, one can construct a cyclic reduction 
preconditioner as in lfl2l . In our tests we used a 3-level V-cycle preconditioner: at the highest 
level, the Schur complement approximation is inverted exactly, whereas for the D rr inversion 
we used vanilla BiCGstab at one percent accuracy. We solved three linear systems to 1CT 8 
accuracy at k = 1/6 on three 8 8 lattices at (3 = 5.4, 5.7, 5.9. The performance of the precondi- 
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Table 1 : Performance of the 3-level V-cycle preconditioner on 8 8 lattices for various couplings. 

tioner is measured using the average convergence rate, 1Z, i.e. the average ratio of the next to 
the current residual error. In the second row of Table Q] we give the estimated condition num- 
bers for each case. In the third and fourth row we show the results for the unpreconditioned 
and preconditioned BiCGstab algorithm. For the first two lattices we used the type I inverse 
permutation, whereas for the third one the type II inverse permutation. It should be noted that 
the preconditioned BiCGstab with the type I inverse permutation did not converge in the third 
case. According to BUl, the lattice at (3 = 5.9 and bare quark mass k = 0.1592 corresponds 
to m n = 320 MeV, which indicates that our pion mass is close to its physical value. This 
shows that, even that close to criticality, a type II inverse permutation will yield a sufficient 
approximation to the Schur complement. 

Applications of this preconditioner in lattice QCD with Wilson fermions are obvious: 

• Acceleration of linear solvers: we didn't make any effort to construct a cheap precon- 
ditioner or to optimise its parameters; our prime interest was to test its potential. A 
production code should take this into account. 

• Acceleration of simulation algorithms: as it is already shown for DD-HMC algorithm 
lfT3~l . the fermion determinant factorisation using the exact Schur complement can effec- 
tively be employed to accelerate the simulation algorithm. A full multigrid preconditioner 
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offers the flexibility to employ more than two factors. 
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